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Abstract. The parameter-free reconstruction of the surface-mass density of clusters of galaxies is one of the 
principal applications of weak gravitational lensing. From the observable ellipticities of images of background 
galaxies, the tidal gravitational field (shear) of the mass distribution is estimated, and the corresponding surface 
mass density is constructed. The noise of the resulting mass map is investigated here, generalizing previous work 
which included mainly the noise due to the intrinsic galaxy ellipticities. Whereas this dominates the noise budget 
if the lens is very weak, other sources of noise become important, or even dominant, for the medium-strong 
lensing regime close to the center of clusters. In particular, shot noise due to a Poisson distribution of galaxy 
images, and increased shot noise owing to the correlation of galaxies in angular position and redshift, can yield 
significantly larger levels of noise than that from the intrinsic ellipticities only. We estimate the contributions 
from these various effects for two widely used smoothing operations, showing that one of them effectively removes 
the Poisson and the correlation noises related to angular positions of galaxies. Noise sources due to the spread in 
redshift of galaxies are still present in the optimized estimator and are shown to be relevant in many cases. We 
show how (even approximate) redshift information can be profitably used to reduce the noise in the mass map. 
The dependence of the various noise terms on the relevant parameters (lens redshift, strength, smoothing length, 
redshift distribution of background galaxies) are explicitly calculated and simple estimates are provided. 
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1. Introduction 

Weak gravitational lensing has been recognized as a pow- 
erful tool to investigate the mass and mass distribu- 
tion ofchKtersof (Webster 1985; Tyson et al. 



199C ; Kochanek 1990 ). In their pioneering paper Kaiser 
fc Squires (1993 ) pointed out that the tidal gravitational 
field ofi tho cluotor, ao moaourod from the diotortion of im 



is that in many clusters the projected mass distribution 
closely follows the distribution of luminous cluster galax- 
ies. 

Recently, a number of papers have been devoted to 
the study of the noise properties of weak lensing mass re- 
constructions. Such studies are important to optimize the 



lensing techniques used. For example, Lombardi & Bertin 



age sh apoe of tho faint background galaxy population, can 
be used to reconstruct the two-dimensional projected mass 
profile s of oluetore, without referring to a paramotoriaod 



1998a) have provided estimates for the noise properties 



of weak lensing mass maps, showing a lso that a particu- 



lar technique, previously described by Seitz fc Schneider 



J 1 Til • ,1 1 1-r; 1 • • 1 ; (1996), leads to optimal maps in the linear regime; anal 
mass n nodol. Ihie method, modmod in various wayf, later — i -j^' i 1: . ti & ' 



(see, e.g.. 


Kaiser 1995; Schneider 1995; 


Schneider & Seitz 


1995; 


Bartelmann et al. 1996|; ^quires & Kaiser 1996; Seit2 


fe Schneider 1997; 3eitz et al. 199S; Lombardi & Bertin 


1998b 


A 


Bartelmann & Schneider 200! 


for a recent re- 



view) has been applied to more than a dozen clusters up 



to now (e.g. pahlman et al. 1994 


Small et al. 1995 


; Squires 


et al. 1996; squires et al. 1996; 


3eitz et al. 1996; 


Luppino 


fe Kaiser 1997; Clowe et al. 1998; Hoekstra et al. 1998; 


MelUer 1999). One of the main 


results from these studies 



ogously, van Waerbeke (2000 ) has extended these results 
to maximum likelihood estimators. 

All these analytical studies, however, have been car- 
ried out using a "mean field" approximation, thus ignor- 
ing some contributions to the noise. In particular, Poisson 
noise and correlations on the positions of source galaxies 
are not taken into account. In this paper we consider in 
detail the effects of these sources of noise and we provide 
analytical estimates for them. We also show that, in many 
situations, these extra contributions to the noise are non- 
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negligible. In partic ular, as already pointed out by Berlin 
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fc Lombardi (2001 ), the noise due to correlations of the 
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source galaxies does not decrease with the density of ob- 
served galaxies, and thus puts an intrinsic limit on the 
accuracy of weak lensing mass maps. On the other hand, 
we also show that the use of some redshift information 
can significantly reduce the noise due to the correlation of 
galaxies and, at least partially, the Poisson noise. 

The paper is organized as follows. We briefly describe 
standard weak lensing mass reconstruction techniques in 
Sect. § In Sect. | we introduce the statistical methods 
used to study the properties of the shear and the mass 
maps obtained in weak lensing reconstructions, and we ap- 
ply those methods to obtain some general results. A very 
useful approximation technique is discussed in Sect. ^. 
This technique allows us to obtain simple analytical re- 
sults and to study the impact of various sources of noise 
on the shear estimation. A more accurate shear estima- 
tor, which leads to significantly smaller noise, is studied 
in Sect. Section |^ is devoted to the impact of obser- 
vations on reconstructions of weak lenses. If information 
on the source redshifts is available, then some sources of 
noise can be significantly reduced. Two different estima- 
tors that can take advantage of redshift information are 
discussed in Sect. ||. Finally, the main results obtained 
in this paper are summarized in Sect. ^ In order to sim- 
plify the discussion and to avoid long calculations, we have 
moved several derivations to the appendices. In particular, 
in Appendix |A we recall a standard statistical technique; 
in Appendices ^ and |c| we carry out the statistical anal- 
ysis for the estimator of Sect. ^; in Appendix ^ we sketch 
the calculations for the estimators of Sect. ^. 

The weak lensing approximation is taken to be valid in 
this paper, so that the shear is taken as an observable [see 
below Eq. (^]. The Einstein-de Sitter cosmological model 
is assumed and the Hubble constant is parameterized as 
Hq = imh km s~^ Mpc~^ 

2. Local shear estimators 

In the weak lensing limit, the observed ellipticity e*^"' of a 
galaxy (n G {1,2,..., TV} denotes the galaxy considered) 
is related to the true source ellipticity e'^t"); 



(1) 



where 7(0) is the lens shear field for fiducial sources at in- 
finite redshift, Z'^ the cosmological weight 
for the n-th galaxy, and is its observed position on 
the sky. We recall that Z{z) depends on the assumed cos- 
mological model, on the redshift of the lens, and on the 
redshift z*^"^ of the galaxy. More precisely, Z{z) can be 
written as a ratio of angular diameter distances between 
the observer and the lens, the obs erver and the source, and 
the lens and the source (see, e.g., Bartclmann & Schneider 



2001) 



average, needed in order to obtain a high signal-to-noise 
estimator from extremely noisy data (the source elliptic- 
ity), is usually performed using a function W{0 — 6') which 
describes the relative contribution of a galaxy located at 
0' to the shear at the position 0. Since only relative values 
of W are important, we can assume that this function is 
normalized, 



/ w{e-e')d^e' = i 



(2) 



where is the field in the sky where the observations are 
made. In principle, the normalization condition written 
above cannot be satisfied for all points 9 oiQ, unless W is 
constant. In practice, if the area of observation 17 is large 
compared to the scale of W , we can assume that Eq. ^ 
holds for nearly all points € 51, being false only for the 
points close to the boundary of the field. This approxi- 
mation, called framing (see Lombardi & Bertin 1998b), 
greatly simplifies the analysis. In the following, we will al- 
ways use this approximation. Often, in weak lensing stud- 
ies, a Gaussian weight function is used: 



w{e-e') = 



exp 



\o-o'\ 



(3) 



Using the smoothing function W, two simple unbiased 
estimators for the shear can be provided. The first one is 
the traditional Kaiser fc Squires (1993 ) estimates for the 
shear: 



m 



1 



N 



:(«)VF(0-6»(")) 



(4) 



Here, p is the density of observed galaxies, for simplicity 
taken to be constant on the whole field of observation (as a 
reference value, a density of 30 galaxies per square arcmin 
can be easily reached today). In this equation and in the 
whole paper we will use the hat {') to denote estimated or 
measured quantities. Equation (Q) is very simple, but has 
the disadvantage of being significantly affected by Poisson 
noise. In fact, it is reasonable to expect some variations of 
the local density of galaxies in the field. In such cases, we 
would have additional noise on 7 simply due to changes 
on p. 

As argued by Seitz fc Schneider (1995| ), since the posi- 



tions of galaxies are readily available, we can modify the 
estimator into 



7(0) 



(5) 



We will refer to this estimator as balanced, because it, 
in contrast to (^), explicitely takes into account the local 
number density of source galaxies. Equation (0) is actually 



often used in weak lensing studies (see, e.g., Lombard: 



Eqi lation ( 1 ) provides a simple way to infer the shear of 
a lens. Since source galaxies can be taken to be randomly 
oriented, so that (e*^) = 0, the average of galaxies close 
on angular position is a good estimator for the shear. The 



ct al. 200C| ) 



In this paper we will study in detail both estimators 
and obtain their expected variances. We will consider ini- 
tially the estimator (M), which is much easier to study than 
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(|). Although non-optimal, this estimator is very impor- 
tant for many reasons: (i) It is the first one considered in 



the literature, in the classical paper by Kaiser & Squires 
(1993|)] (ii) it is simpler to study and hence allows exact 
analytical results; (iii) it has been used in the literature 



(e.g. [Fahlman et al. 1994 


Tyson & Fischer 1995) and ac- 


tually is still used by some authors (Fischer et al. 1997; 


Fischer & Tyson 1997; I 


juppino & Kaiser 1997 




Fischer 



hoping also that a detailed analysis of the unbalanced es- 
timator will make the complex discussion of the balanced 
one easier to follow and results interpreted. 

An estimate of the shear can be used to obtain an es- 
timate for the projected, dimensionless mass distribution 



& Schneider 1996: Lombardi 



e.g., 


Seitz 


1999a 


). In 



general, because of the linear relationship between 7 and 
K, we can write 



K{e)^K + / A(0,e')7.(e')d 



2/)/ 



(6) 



where the Einstein's convention on repeated indeces has 
been used. In this equation, R is the mean value of k, on 
ri and 'Di{9,6') is a suitable kernel. For example, in the 
limit where the area of observation is much larger than 
the size of the lensing cluster, a simple kernel can be used 
(Kaiser & Squires 1993): 



Note that in this case Eq. (|^) is a simple convolution. For 
finite fields fi, the kernel l?i(0, 9') often cannot be written 
in closed form, still it can be written in terms of the Green 



functio n for a Neumann proble m (see Lombardi & Berlin 
1998a[n $eitz' fc Schneider 200l| ). Fast numerical methods 
exist to obtain k from the shear 7 ( Lombardi fc Berlin 
1999a|)l 

The linear relationship between k and 7 allows us to 
easily derive the statistical properties of the mass map 
from the properties of the shear. In the following, thus, 
we will first evaluate the expectation value and the noise 
of the shear, and then obtain analogous results for the 
mass map. 

3. Statistical analysis 

In this section we will perform a statistical analysis of 
the shear estimator (^ using a general probability density 
distribution for all random variables involved. As already 
noted above, this estimator is not optimal and generally 
should not be used. One purpose of this paper is actually 
to show that the increase in noise of unbalanced estimator 
is very large in most cases. 

Let us describe in detail the various random variables 
used. 

Source ellipticities e''("). We assume that the intrinsic 
ellipticity of each galaxy follows the probability distri- 
bution pg{e^); this probability distribution is taken to 



be isotropic, so that (e''(")) = 0. Moreover, we assume 
no correlation between ellipticities of different galax- 
ies, and no correlation between ellipticities and other 
random variables. 

Source redshifts 2*^"^ The redshifts of sources are 
taken to be unknown in the first part of this paper (see 
Sect. ||for the case where some knowledge on redshifts 
is available). We assume that each redshift follows the 
probability distribution ^^(z). Correlations in redshifts 
and between redshifts and positions are considered. 

Observed position S^^"^. The positions of galaxies are 
obviously readily available in weak lensing observa- 
tions. However, in order to obtain general results, we 
will perform an ensemble averaging on galaxy po- 
sitions. This way, we will obtain results which are 
not strictly dependent on the particular configuration 
of galaxies. Here we will assume that observed posi- 
tions follow a uniform distribution on the field of ob- 
servation (thus neglecting the magnification bias). 
Correlations in positions and between positions and 
redshifts are considered. 

In summary, the average of a function / of random vari- 
ables of a single galaxy will be performed using the relation 

(/(e("),z("),0("))) = l|^d20^1(z)dz 

X / p,{e')f{e' + Z{zh{e),z,e)d'e^ . (8) 

|e=t<l 

Here A is the area of the observation set fi. Similarly, 
the average of the product of two functions / and g of 
random variables of two different galaxies (n 7^ m) will be 
performed using the relation 



/(e("\z("),6>("))\ • (.g(e(™\z("'),6>("')) 



+ 72 / dV / p.{z)dz / p,{z')dz' 

^ Jn Jn Jo Jo 



X / p,ie^)dh^ f p,{enfie,z,e)g{e',z\e')^d\^' 

J|e''|<l J|e="|<l 



(9) 



In this relation, for simplicity, we have kept the original 
arguments e = e'' -I- Z{z)-f{e) and e' = e"' + Z{z')j{e') for 
/ and g in the last term. The function ^, for our purposes, 
can be written as 



^ = a0-9\z,z') 



(10) 



The galaxy two-point correlation function is often written 
in terms of the proper distance d between the galaxies. 
Observations show that the local galaxy two-point cor- 
relation function is very well approximated in the range 
10 kpc < hd < 10 Mpc by the power law 



(11) 
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with do — (5.4± 1) h^^ Mpc, and rj ~ 1.77±0.04 (Peebles Analogously, we define the two-point correlation function 
1993). At high redshift (z ^ 1), the two-point correlation for the mass map as 
function is more poorly known. The redshift-dependence 
is often written as a power-law 



(12) 



Recent investigations ( Le Fevre et al. 1996 ) suggests that 
a is in the range [—3, —1]. In this paper, then, we will take 
the fiducial value a = —2. 

Finally, we note that, given the linear relationship be- 
tween e and e^, it is not difficult to include the effect of 
errors of measurements for the observed ellipticities in our 
discussion. In general, such errors can be described using a 
Baysian probability distribution Pobs(e|e) which gives the 
probability of measuring e when the real lensed ellipticity 
is e. Then the standard framework described above can 
be applied without any change if we replace the probabil- 
ity distribution p(e*^) with an effective distribution poff(e^) 
given by 



Poff(e') 



Pobs(e1e^)p.(e^)d2e^ 



(13) 



e= <1 



We stress that, because of the linearity of Eq. (|l|), the dis- 
tribution Pobs(e^|e'') is formally identical to the distribu- 
tion Pobs(e|e)- Finally, we observe that if this Baysian dis- 
tribution is taken to depend only on the difference |e — e|, 
then the effective variance cr'^^g of e*^ can be written as 

^eoff = f^e + <^lhs > (14) 

where a^^^ is the variance of (e — e) (given by Pobs)i and 
(T^ is the variance of the source ellipticities, defined as 



is easily evaluated 



3.1. Average value 

The average value of the estimator 
using Eq. (^): 

(7(0))= / wie-e')-f{e')d^9' . (15) 

Jn 

Because of the linearity between 7 and k we also obtain 

{k{e))= f wie-e')Kie')d^9' . (i6) 



In other words, the expected value of the estimated mass 
map is the true mass map smoothed by the weight func- 
tion W. The length-scale of this function, thus, sets the 
resolution of the shear and mass map. 

3.2. Covariance 

A complete characterization of the error on the shear map 
is given by the two-point correlation function, defined as 
(see Lombardi fc Berlin 1998a) 



Gov,, (7; 0,0') 



[7.(0) -(74^))] [ij{d')~{ij{e'))] 



(17) 



Equation 
Cov(7): 

Cov(k;0,0') 



(18) 

provides a simple way to obtain Cov(k) from 



d20i?,(0,</,) / d^^'v,{e\ci,') 



X Covy(7;0,0') 



(19) 



We stress that the simple variance of 7, i.e. the quantity 
Covii(7, 6, 9), would not be sufficient to evaluate the vari- 
ance of k. 

In order to obtain Cov(7) we write the scatter of 7 
from its average value as 

1 ^ r 

[7. - (7.)] (^) = 7^ E (^^^"^ + 1^{0^''^)Z^''^)W{9 0(")) 



n=l 

{Z)_ 
A 



w{e-ci>)-f,{ct>)d^ 



(20) 



since N = Ap. Hence we find 
Gov,, (7; 0,0') = 

/ AT - 

1 



(Z)2p2 



^ (ef^) -f 7,(0("))Z("))l^(0-0(") 

{Z) 



Y / W^(0- 0)7.(0) d'' 



N p 



s(rn) 



^ |^M/(0'-0')7,(0')d2 



(21) 



This expression is evaluated more easily if we distinguish 
the two situations n = m and n 7^ m. Let us start with 
the case n — ni. 

If n = m, then we are considering the same galaxy and 
Eq. (j|) must be applied. In this case, using the isotropy 
hypothesis for source ellipticities, we obtain 



{Z?P 

1 



/ iy(0 - 0)iy(0' - 0) d2(/- 

7.(0)7,(0)m^-0)VF(e' 

I 

M^(0- 0)7,(0) d2</) 
W^(0'-0')7,(0')d20' . 



0)d^ 



(22) 



The last term vanishes in the limit where the area of ob- 
servation A is large. [What actually matters here is the 
ratio between the "effective" area of the weight function 
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W and A.] We thus will ignore this term and similar terms 
arising in other expressions. 

When n 7^ TO, on the other hand, we obtain a single 
term involving the two-point correlation function 



^1 Jn Jn 



C. = 



X / p,{z)Z{z)Az / p,{z')Z[z')^{(l)-(t)',z,z')<lz' . 
Jo Jo 

(23) 

The last two integrals in the r.h.s. of this equation are rem- 
iniscent of Limber's equation (see Peebles 1993| ) and could 
suggest that the (projected) angular correlation function 
is relevant for lensing. This function, in fact, can be writ- 
ten in term of the full three-dimensional correlation ^ as 



Pz{z)dz / pz{z')£_{(p,z,z')dz' 



(24) 



In reality, because of the presence of the cosmological 
weights Z{z) and Z{z') in the integrand of Eq. ([2^), the 
expression for C2 cannot be rewritten in terms of the an- 
gular correlation w(ip) only. 

To simplify the notation, we will write terms similar 
to the one found in Eq. (03) as 



C2 = 



1 



{zy 



■E[W-,,ZW'j'jZ'] 



(25) 



The functional S represents a multiple integration over 
^{(j)—ip',z,z') d^0d^0' dz dz' . In summary, the correlation 
for the shear can be written as 

Gov,, (7; e, e') = j0- w{e - ct,)w{e' - 0) d^^ 



1 



/ -i,{ct>)^,{ct>)w{e - 4>)w{e' ~ ct>)d^ 

Jn 



E\^-i^ZW'i^Z'\ 



(26) 



This expression, although rather complicated, already 
clearly shows the various contribution to the noise. The 
first term on the r.h.s., proportional to g^, is due to in- 



trinsic scatter of source ellipticities. Lombardi & Bertin 
(1998a|) have already considered this term, and have also 
shown that because of its "diagonal" form (it is propor- 
tional to 5ij and is a simple convolution of W with itself) 
can be easily converted into a covariance for k. Writing 
(Z^) = {ZY + cr^, the second term on the r.h.s. can be 
thought as the sum of two terms. One term does not in- 
volve any redshift dependence and is a Poisson noise on the 
angular density of galaxies. The other term is proportional 
to the scatter of galaxies in redshift, and is a Poisson 
noise on the redshift distribution. These two terms can 
be better understood by two simple examples. Suppose 
that we observe some gravitational lens and that we make 
a shear map using the estimator of Eq. (^). Because of 
Poisson noise, in some region of the sky we could observe 
an overdensity of background galaxies, and thus obtain a 



larger than expected value for the shear in that region. 
On the other hand, in another region we could observe a 
"normal" density of galaxies, but their redshifts, because 
of Poisson noise, could be larger than expected. As a re- 
sult, in this second region too we would over-estimate the 
shear. Situations such as the ones described contribute to 
the noise on 7 by adding two sources of errors character- 
ized by their dependence on 7 (in both cases considered 
above, in fact, no error is expected if 7 = 0). Finally, 
turning to the last term of Eq. (p6[), we note that it is 
clearly related to the two-point galaxy correlation func- 
tion. Analogously to the second term of Eq. ([2^), this is 
actually composed of two contribution, one due to the an- 
gular correlation and one due to correlation in redshift [cf. 
Eqs. (H) and (H)]. 

From Eq. ( p6| ) we can also see that the first two noise 
terms are proportional to while the last one is appar- 
ently independent of the density of background galaxies. 
This behavior is easily explained. EUipticity and Poisson 
variances depends on single galaxies and thus increase lin- 
early with p; hence, the factor 1/p in Eq. (^) makes the 
final noise due to these sources proportional to 1/ On 
the other hand, the variance due to galaxy correlation is 
proportional to the number of pairs of galaxies and thus 
to the final correlation noise is then formally indepen- 
dent of p. We stress, however, that all noise terms depend 
on the galaxy redshift distribution. The density p and the 
probability distribution Pz{z) are intimately related. For 
example, deeper observations will increase p but also the 
mean galaxy redshift. As a result, all noise terms do de- 
pend on the depth of the observation (see Sect. ^ for fur- 
ther comments on this point). 

Looking again at Eq. (p6|), we can also obtain informa- 
tion about the typical correlation length for the shear, i.e. 
the largest angular separation between points that have 
correlated noise properties. The first term in this expres- 
sion, in fact, suggests that the correlation length for 7 is 
approximately twice the length scale of the weight func- 
tion W; the same is true for the second term, while the 
third term has no "intrinsic" scale. This last point is basi- 
cally due to the lack of a scale of the two-point correlation 



function as given by Eq. (11) (note that do in this expres- 
sion should be taken as a simple multiplicative coefficient 
rather than as a scale for ^). 

In principle, for a given lens, we could now evaluate the 
covariance of 7 using the results obtained so far, and then 
convert it to the covariance of k. Actually, in the limit of 
large fields fl, we can use the kernel (^ in Eq. (^9|), thus 
obtaining directly an expression for the covariance of n: 



Cov(k; 6>, 6»') 



{zyp 



w,ie ~ c(,)w,{e' ~ cf>) d^ 



{Z)^P Jn 
1 



{zy 



7,(0)7,(0)1^,(6/ - <i))Wj{e' - 0) d^(j) 

E[W,-i,ZW'ji^Z'] . (27) 
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Here, Wi is a "modified weight," i.e. the convohition of W 
with X>^S. 



0)i?P(0)d^0 



1 



exp 



(28) 



The last equality holds when a Gaussian weight function 
of the form (0) is used. 



4. Slowly varying shear 

The results obtained in the last section are exact but, un- 
fortunately, difficult to interpret in general. In particular, 
it is not easy to estimate the order of magnitude of the 
expressions obtained, and it is also difficult to say which 
term makes the more important contribution in Eq. ( ^6|) 
[or, equivalently, Eq. (27)]. In an important limit, how- 
ever, the results obtained so far can be written in a much 
simpler and more useful form. 

Suppose that the true shear field 7(6) is slowly vary- 
ing on the typical scale of the weight function W . In this 
limit, we can expand the shear to first-order and greatly 
simplify the integrals that involve 7 and W . The results 
obtained, though clearly not exact, represent an extremely 
useful starting point to evaluate the noise of the proposed 
shear estimates. Note that this technique can be success- 
fully used to simplify the expressions of statistical esti- 
mates for the shear field, while is not directly applicable 
to expressions for the mass map because of the non-local 
dependence of k on 7 (or, equivalently, because of the 
infinite support of V^^). Actually, we will still be able 
to obtain a finite form for the average value of k{6) (see 
Sect. but not for the covariance. However, if the noise 
properties of the shear are known, the covariance of k can 
be evaluated using Eq. (|l9|). Moreover, it is not unreason- 
able to assume that the covariance of k is close to Cov(7) 
[actually, if only the first term of Eq. ( p6| ) is taken into 
account, the two covariances are identical]. 

In the following we will always assume a symmetric 
weight function, i.e. such that W{0 — 6') depends only 
on \0 — 6'\. Moreover, the results obtained will also be 
specialized to the important case where is a Gaussian of 
the form (^) . Finally, in order to numerically estimate the 
contribution of some terms to the noise, we will often refer 
to a "typical" model. The model is characterized by a lens 
at redshift = 0.3 in an Einstein-de Sitter cosmological 
model. The normalized redshift probabili ty distribution 
of so urce galaxies will be taken to be (see Brainerd et al 



19961 ) 



Pz{z) 



exp 



The constant zq determines the depth of the observations, 
while P determines the steepness of the cutoff at high red- 



s' 



Fig. 1. Linear terms of the expansion of 7 around 9 do 
not enter the expression for the covariance. In fact, the 
integrand of any linear term evaluated at the point (pi is 
the opposite of the integrand evaluated at (p2 {<i>i and 4>2 
are such that (pi + (p2 = 20). In particular, the applica- 
tion of this symmetry to the liner expansion for the term 
7,(0)7j (^) leads to Eq. (||). 



shifts. The moments of this distribution are 



and the mode is 



■^modc — ^0 



1//3 



(31) 



In the following, when not otherwise stated, we will use 
P — 1.5 and zq = 0.7, so that (z) ~ 1.05 and Zmodc — 
0.848. 

4.1. Average value 

If 7(0') does not change significantly when 9' is on the 
effective area of W, we can expand the shear around the 
point 9: 



7.(0')^ 7(0) 



9,) 



(32) 



Inserting this expression into Eq. ([Tsl), and using the sym- 
metry oiW, we get 



im) 



7(0) / W{9 
Jn 



9')d^e' ^-l{9) 



(33) 



Equation ( |33|) is clearly a rough approximation, since 
we already know that average value for the measured shear 
map is a smoothed version of the true shear. On the other 
hand, if the shear does not change significantly on the scale 
of W , we can safely take the result obtained; in more gen- 
eral cases it is, at least, an order-of-magnitude estimate. 



(29) 4.2. Covariance = 0) 



We can in a similar way obtain an expression for the co- 
variance of the shear using Eq. (|2^). The method used 
is similar to the one adopted above, with the only small 
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difference that now the expansion for 7 will be performed 
around the mid-point 6 = {e + e')/2 (see Fig. |l|). We thus 
find 



Cov,j(7;0,0') 



+ 



X (w*w){e^e') 



lr,A0)l,.3'{e) 



X / w{e-<t,)w{e' ~<t,){(t),. 

In 



(34) 



where the star (★) represents a convolution. In the case 
where is a Gaussian function we obtain in particular 



Gov,, (7; 0,0') = 



1 



{Z-)a 



X exp 



\e-e 



,/|2 



(35) 



This is the result for the covariance of 7 in the slowly vary- 
ing shear approximation. The fact that a finite expression 
has been obtained allows us to study the behavior of the 
various contributions to the noise. Let us discuss in detail 
the result obtained. 

Equation (|3j) shows that the noise is proportional to a 
convolution of W with itself. This suggests to define the ef- 
fective area A\y of the weight function, i.e. the area where 
W{(f)) is significantly different from zero or, formally. 



A 



w 



(36) 



For the Gaussian weight function we get Aw — 47rcr^ 



We then note that the noise on 7 given by Eq. (|3^) is 
proportional to 1/iVw, where Nw = pAw is the expected 
number of galaxies inside the effective area of W . 

The ratio between the contributions to the noise aris- 
ing from the shot-noise and the intrinsic scatter in ellip- 
ticities is of order \^{0)\ (Z^)/cr^. Since — 0.25, and 
(Z^) ~ 1, we see that the shot-noise is a relevant term 
unless the lens is very weak. More specifically, let us con- 
sider the error made on the component i of the shear at 
6. We have from Eq. (p5h (no sum on i) 



Gov,, (7; 0,0) 



1 



al + {Z^)[^mY 



^po ■ 



y2 

^Pl ; 



(37) 



where, for convenience, we have explicitely split the three 
contributions to the noise considered here. Using the lens 
model described above, we have {Z) ~ 0.600810 and a% ~ 



0.038347. Hence, the expected error on 7 due to the spread 
of source ellipticities is 



= 0.016600 



0.25 7 V 1 arcmin / 

-1/2 



50 gal arcmin 



(38) 



The Poisson noise contributes with a term proportional to 
the shear, 

Spo = 0.004196 (2_Y^^Y 

^ rm / 



0.1 / V 1 arcmin^ 

1/2 



50 gal arcmin 



(39) 



and with a term proportional to the derivatives of the 
shear, 



Spi 0.002967 



7. 



0.1 arcmin / \50 gal arcmin 



P 



-1/2 



(40) 



Note that, surprisingly, this term does not depend on aw 
(in this case, the angular scaling is given by the gradient 
of the shear). In conclusion, for a typical weak lens at 
redshift Zd — 0.3 we expect similar contributions from the 
intrinsic scatter of ellipticities and Poisson noise. 

4.3. Covariance (S, 0) 

Analytical calculations for the contributions to the covari- 
ance from the galaxy two-point correlation are definitely 
non-trivial. For this reason, we report here only the re- 
sults of such calculations, and we refer to Appendix |^ for 
any details. Note that the calculations have been carried 
out only in the (important) case where is a Gaussian 
of the form (^). Moreover, sensible approximations have 
been used. 

As discussed above [see Eqs. (|l^) and (^], the two- 
point correlation function is well approximated by a simple 
power-law in the physical separation of galaxies. Moreover, 
observations suggests that the exponent 77 is very close to 
2. This allows us to carry out calculations, obtaining [see 
Eq. (H)] 



Co 



exp 



{i + zr[p,{z)z{z)y 

Da{z)D'p(z) 



Az 



(41) 



where /q is the modified Bessel function of first kind, ip = 
6' 6, and Da{z) and Dp{z) are, respectively, the angular 
diameter distance and the proper distance of an object at 
redshift z [note that in Eq. (41) the derivative of Dp{z) is 
involved]. Note that, in Eq. (41), the smoothing length aw 
has to be expressed in radians. Moreover, we have used as 
usual the notation for the mid-point between and 0'. 
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In the typical case discussed here of a lens at redshift 
Zd = 0.3 and using do — 5.4 /i"^ Mpc we have for the 
expected error due to the correlation of galaxies, S2 = 



E2 = 0.031173 



_L 
0.1 



\ -1/2 

1 arcmin / 



(42) 



We must thus deduce that, at least for the case considered 
here, the galaxy two-point correlation function represents 
a major source of error in weak lensing estimates. A more 
detailed discussion on the relative contributions of the var- 
ious sources of errors is delayed until Sect. |^. 

5. The balanced shear estimator 

As already pointed out, the estimator is very simple 



and traditionally has been the first one used (Kaiser & 
Squires 1993). On the other hand, this estimator has a 



rather large noise and in particular suffers from Poisson 
noise, due to the fact that the normalization factor is con- 
stant both on the density of galaxies (1/p) and on the 
redshift distribution of sources (1/(Z)). 

In this section we will study in detail the estimator (ph 
similarly to what we have already done in Sects. |3| and^ 
for the estimator (^). However, since now the estimator is 
non-linear on the observed positions of galaxies, an ana- 
lytical evaluation of the average values and errors for this 
estimator would be hopeless without any approximations. 
For this reason we will linearize Eq. (|^) around the av- 
erage values of the random variables involved. This tech- 
nique, often adopted in statistics, is briefly described in 
Appendix^. Moreover, to simplify the discussion, we will 
mostly report here only results, referring to Appendices ^ 
and |c| for the derivations. 

5.1. Linear approximation 

Let us call N and D the numerator and denominator of 
Eq. (||), so that 7^ = Ni/D. Using the linear approxima- 
tion (see Appendix]^, we can write the expectation value 
of 7i as 



(7. 



The last two terms on the r.h.s. represent second-order 
corrections and will often be ignored. The covariance of 
7 can be written as Cov^ (7; 9, 9'), i.e. as a 2 x 2 matrix 
which depends on two spatial coordinates 6 and 6' (see 
Lombardi & Berlin 1998a). Omitting the two arguments 



9 and 9' , we can write 



Cov,y(7) 



1 



{D){D'] 



■ CoviN^N'] 



{D){D') 



■ Cov{N,D') 



{Dy{D' 



■Cov{DN') 



{D)^D'Y 



Cov{DD') 



(44) 



Here we have used the prime to denote quantities which 
must be evaluated at 6' rather then at 6. 

As a general rule, the linear approximation can be used 
when the denominator D is expected to have a small (rel- 
ative) variance. This is true if the expected number of 
galaxies inside the effective area of the smoothing func- 
tion W is much larger than 1. 

We also note that in the linear approximation, the ex- 
pression for the covariance for 7 depends on the two-point 
galaxy correlation function only, and not on higher-order 
correlation functions. This would not be true without the 
linear approximation. 

5.2. Slowly varying shear 

Similarly to what we have seen above for the unbalanced 
estimator (^ , the statistical properties of the balanced es- 
timator (|^) are better understood in the limit where the 
shear is slowly changing on the effective area of the weight 
function. In this shown in detail in Appendix |C| 

we can explicitely write the expectation value and the co- 
variance of 7. 

Regarding the expectation value of 7, we obtain a re- 
sult similar to Eq. (p^), i.e. (7(0)) — ^{9). In other words, 
in the slowly varying shear approximation, the estimator 
(^ is unbiased. Note that in more general cases we still 
recover Eq. (|l5|), i.e. the average value for the estimator 
7 is the true shear 7 convolved with the weight function^. 
Finally, since Eq. ( |l5[) holds, we can immediately deduce 
that Eq. (16) holds as well. 

Turning to covariances, we can finally see the advan- 
tages of using a balanced estimator for the shear. In fact, 
in the slowly varying shear approximation and neglecting 
the two-point correlation function, we obtain for the co- 
variance of 7, in the case where is a Gaussian function. 



Cov,j(7;0,0') = 



w 



X exp 



2 

Ml 
4-^ 



(45) 



with cr| = (Z^) - (Z) . This expression should be com- 
pared with the analogous equation obtained for the esti- 
mator (|^), namely Eq. (35). The basic difference is given 
by the factor (Z"^) which replaces cr| of Eq. ( |4^ ) or, 
in other words, by the presence of an extra term (Z)'^ 
in Eq. (|35|). This extra term is precisely due to Poisson 
noise on the angular coordinate 9; analogously, the term 
cr| left can be interpreted as Poisson noise in redshift. 
In order to better appreciate this point, it is interesting 



^ Actually, a rigorous calculation shows that the average 
value for 7 is the true shear 7 convolved with a kernel slightly 
different from W . The kernel quickly converges to the weight 
function when the number of effective galaxies is large (see 
[Lombardi fc Schneider 200 1[ ). 
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to evaluate the expected error of the balanced estimator 
in the lens model described above. In this case we have 
{Z) ~ 0.600810 and cr| ~ 0.038347. Hence, the use of a 
balanced estimator reduces the "constant" component of 
Poisson noise by about a factor 4 (i.e. a factor 16 in vari- 
ance). More explicitely, the contribution of the redshift 
Poisson noise is given by 



Spo = 0.001300 ( ^ 



P 



1 arcmin / 

-1/2 



50 gal arcmin 



(46) 



This expression replaces (|39|). Note that, in contrast, the 
Poisson noise due to shear variations [represented by the 
second line in Eq. (^5|)], is only partially reduced by the 
balanced estimator. In particular, if y = 0, no improve- 
ment is obtained and Eq. (^0|) can still be used to have an 
estimate of this noise source. 

If the two-point galaxy correlation is taken into ac- 
count using the model (^2|) , we have a supplementary term 
in the covariance. For the lens considered here, this term 
is of the order (see Appendix |c|) 



0.007152 ( ^ 



1 arcmin / 



-1/2 



(47) 



Comparing this result with Eq. ([42[), we see that the use of 
the balanced estimator has significantly reduced the noise 
due to the clustering of galaxies. 

6. Numerical estimates 

In the previous sections, we have obtained estimates for 
the various contributions to the noise of the shear using 
several approximations, namely the weak lensing limit, 
the slowly varying shear (Sec t. |4| ), and the linear expan- 
sion of the estimators (Sect. |j). Among these approx- 
imations, probably the most critical one is that of the 
slowly varying shear. Hence, in order to test the relia- 
bility of the results obtained so far, we have numerically 
evaluated the noise on the shear estimator (||) in a spe- 
cific case. Note that, instead, it is extremely difficult to 
test directly the results for the balanced estimator with- 
out using the linear approximation described in Sect. 5 
In this case, in fact, we would need to average over the 
whole galaxy positions using also a probabihty 

distribution p[9'^^\ 6^'^\ . . . , 6^^^) which includes effects 
of n-point correlation (actually, even TV, the total number 
of galaxies, should be taken as a random variable with 
Poisson distribution). Clearly, such an analysis would be 
prohibitive. 

We have considered a lens at redshift — 0.3 with an 
axisymmetric profile (see Fig. |2]): 



<0) 



-f{\o\/o. 



-f{\O\/0t 



(48) 



where f{x) is the dimensionless function 



(49) 



Note that f{x) (x l/x for large x, k.{6) cx l/\9\ for 0^ < 
\e\ < Ot, while K{e) cx l/\e\^ for \e\ > Ot. For this lens, 
the shear is easily obtained: Its modulus is 



l{0) 



with 



-gim/Oc 



KqOc 



-9{\e\/0t 



(50) 



(51) 

for large \6\ (see 



The shear decreases as j{6) oc 1 
Fig. I). 

Using this lens model, we have carried out the inte- 
grations of the various terms of Eq. ( p6| ) numerically, and 
compared the results obtained with the analytical predic- 
tions of Sect. H Table ^ shows the relevant parameters 
used in the calculations. The results obtained are summa- 
rized in Figs. ^ and ^. 

Figure || shows the lens profile together with the tan- 
gential shear. In this figure we have also plotted the weight 
function W{6) used in the calculations. As shown by 
Fig. 1^, the slowly varying shear approximation is not really 
justified in this case, since the lens shear shows significant 
deviations from linearity on the scale of the weight func- 
tion. However, we find interesting to consider this situa- 
tion in order to better test our analytical approximations. 

Figure ^ shows, for different radii, the expected er- 
rors on the shear coming from different sources of noise. 
The horizontal line is the intrinsic cUipticity noise, which 
clearly is constant on the whole field. Note that the ana- 
lytical estimate for this noise source is exact. As shown by 
Fig. 1^, correlation of galaxies is the main source of noise at 
\0\ ~ 1', while Poisson noise is relatively unimportant (see 
Sect. below for further comments on this point). The an- 
alytical approximations given in Sect. ^ perform surpris- 
ingly well, with deviations below 10% for most radii. The 

Table 1. Relevant parameters used for the numerical in- 



Comment 



tegration. 




Variable 


Value 


Zd 


0.3 


Oc 


30" 


et 


5' 




0.8 


P 


50 arcmin" 




0.25 


zo 


0.7 


P 


1.5 


do 


5.4/1 Mpc 


V 


1.77 


a 


2 


aw 


20" 



Lens redshift 

Lens core radius [Eq. (BS) 
Lens truncation radius [Ec 
Lens central density [Eq. (4 
Density of galaxies 
Intrinsic ellipticity scatter 
Redshift distr. parameter [Eq. 
Redshift distr. slope [Eq. (|29|)] 
Galaxy correlation length [Eq. 
Galaxy correlation slope [Eq. ( 
Galaxy correlation redshift 
dependence [Eq. ([l^)] 
Smoothing length 
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Fig. 2. The lens mass profile k{6) and the lens shear 7(0) 
at difi^erent radii. The plot also shows the Gaussian weight 
function W{0) used. Notice that w{9) in this plot is not 
normalized. 
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Fig. 3. The different contributions to the shear error (i.e., 
square root of variance) at different radii. Note that the 
analytical estimates, obtained from the expressions of 
Sect. ^ are in very good agreement with the numerical 
results. 



" 0.06 



n 1 

ellipticity 
Poisson num. 
Poisson an. 
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Fig. 4. The different contributions to the mass noise at 
different radii. Although the analytical curves refer to the 
shear noise, they provide a very good estimate for all radii 
larger than 1'. The scale of this figure is the same as Fig. |[ 
the curve for the numerical evaluation of the correlation 
noise reaches about 0.67 at the origin. 




Fig. 5. The mean cosmological weight and its moments as 
a function of the lens redshift z^. Note that the quantity 
(7z fia-s a peak around Zd = 0.5. 



larger deviations observed around \6\ — 1' are expected, 
since at this radius the shear deviates significantly from 
the linear approximation (the same is true around the ori- 
gin, where the shear vanishes). We also point out that the 
analytical estimate for the correlation noise is slightly less 
accurate than the one for Poisson noise mainly because 
Eq. ( ^ ) is a zero- order estimate (in the sense that terms 
containing shear derivatives have not been considered in 
this equation). 

Figure |^ shows similar results for the noise on the 
mass map k{6). The numerical integrations have been per- 
formed using Eq. (27), while the analytical estimates are 
the same as for the shear (i.e., we have assumed that the 
noise on the mass map is the same as the noise on the 
shear). Note that, except at small radii, the analytical es- 
timates for the shear provide an accurate measure of the 
noise on the mass map. In reality, for more complicated 
lenses this property is not verified at the degree of pre- 



cision of Fig. U; nevertheless, the analytical expressions 
for the shear noise are always a good order-of-magnitude 
estimate for the noise in the mass map. 

7. Impact on observations 

So far we have studied the statistical properties of two 
shear estimators from a purely analytical point of view . 
Clearly, the interest of our results lies on the consequences 
of the various noise sources on weak lensing studies. In this 
section we consider the expressions obtained in Sects. ^ 
and ^ from a more practical point of view, with particu- 
lar attention to their dependence on critical observational 
parameters (such as the redshift and the strength of the 
lens, the density of galaxies, and their mean redshift). 

As already stressed, various sources of noise contribute 
to the expected error of both estimators. In particular, we 
have found three main contributions: 
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Fig. 6. Typical errors on the shear due to various sources 
of noise. The plot shows the numerical factors in front 
of some equations describing errors on 7. For all graphs 
we assumed p = 50 gal arcmin"^, aw ~ 1 arcmin, and 
7 = 0.1. For the estimator (Q), Eq. ( p8|) gives the error 
due to the intrinsic scatter of source cUipticities (as- 
suming (Tj = 0.25), Eq. the Poisson noise Epo (both 
in redshift and in position), and Eq. (42) the noise due to 
clustering of galaxies S2; for the estimator (^), the noise 
due to source elhpticities is still given by Eq. (^, the 
Poisson noise (now only on redshift) is given by Eq. (^6|), 
and the clustering noise by Eq. (p7|). Note that the relative 
importance of error source depends critically on the shear 
7. For example, for a lens with 7 = 0.3 the clustering noise 
for the balanced estimator would be larger than the intrin- 
sic ellipticity noise in the redshift range 0.2 < Zd < 0.9. 



1. the intrinsic scatter of source elhpticities; 

2. the Poisson noise, both on the positions of galaxies and 
on their redshift; 

3. the clustering of galaxies, characterized by the two- 
point correlation function ^. 

Depending on some key parameters, these sources of noise 
can have different relative importance for the total error 
of the measured shear, and thus of the reconstructed lens 
mass map. 

For weak lenses the leading term is the intrinsic scatter 
of elhpticities. In fact, the Poisson noise and the noise due 
to the clustering of galaxies are proportional to the shear, 
and thus are expected to have a negligible impact on very 
weak lenses. We note that this point is very important for 
cosmic shear studies. 

The noise due to galaxy clustering shows rather sur- 
prising properties if compared to other sources of errors. 
As already pointed out in Sect. 3.2, ellipticity and Poisson 



noise decrease with the density of background sources as 
contrast, the noise due to clustering is formally 
independent of the density of galaxies. Hence, this source 
of error is expected to be particularly important in deep 
observations and can represent a major limitation for weak 
lensing studies in the near future. In reality, deeper obser- 
vations also imply a change of the redshift distribution 
and, in particular, an increase of (z). As a result, the cor- 



relation error decreases for deep observations, but not as 
fast as other sources of noise. 

Table ^ reports typical figures for the various noise 
terms as a function of the density of background galaxies. 
Calculations for that table have been carried out using 
(Te = 0.25, aw — 1 arcmin, and a lens shear 7 — 0.1. 
In order to simulate the increase of the mean redshift for 
deep observations, we have used a galaxy redshift prob- 
ability distribution Pz(z) of the form of Eq. (29), with 
f3 ~ 1.5 and zq dependent on the density p. In particu- 
lar, we have used zq = 0.6 for p = 25 arcmin^^, zq = 0.7 
for p = 50 arcmin^^, and zq = 0.8 for p = 100 arcmin"^. 
Looking at Table |^, we find surprisingly that the ellipticity 
noise decreases as p-O-SO instead of as p~^/^, as suggested 
by Eq. (^8|). The 0.1 increase in the logarithmic slope of 
Se is due to the change in the galaxy redshift distribution 
considered here, and in particular is related to the coef- 
ficient 1/{Z)^ in Eq. (|3^). In contrast, the Poisson noise 
Epo for the unbalanced estimator almost exactly keeps 
the p~^l'^ dependence. This is due to a cancellation effect 
in Eq. (|37|), between (Z)^ in the denominator and (Z^) 
in the numerator. The correlation noise S2 for the unbal- 
anced estimator is found to be basically constant with p; 
again, the 0.1 difference between the expected logarithmic 
slope (S2 c>c p°) and the observed one (S2 oc p^°'^) can 



be attribute to a 1/(Z)^ factor [see Eq. (41)]. 

Turning to the balanced estimator, we immediately 
note a more pronounced dependence of the noise terms on 
p. The Poisson noise decreases as p~°'^^, i.e. with a loga- 
rithmic slope 0.22 larger than suggested by Eq. (|4^). We 
can certainly attribute 0.10 increase to the factor 1/(Z)^ 
in Eq. (^5|); the remaining 0.12 increase is due to a\., which 
decreases as the bulk of sources is shifted towards high 
redshift. The correlation noise is found to have a density 
dependence p~^-^^ . It is interesting to note that for both 
Poisson noise and correlation noise we gain an exponent 
—0.2 with respect to the unbalanced estimator. 

For both estimators, the correlation noise is found 
to be significant (we recall that calculations for Table |^ 
have been carried out assuming a shear of only 0.1). 
Fortunately, as we will see in the next section, this source 

Table 2. Noise terms as a function of the observation 
depth. The table reports the expected noise contribu- 
tions (square root of variance contributions) from various 
sources of noise. The last column gives the logarithmic 
slope of the respective noise term in terms of p, i.e. the 
constant fc in the relation Ex oc p^ . 



Noise term 


Density p (arcmin 
25 50 100 


k 


Ellipticity noise Ee 


0.025689 0.016600 
Unbalanced estimator 


0.010996 


-0.60 


Poisson noise Epo 
Correlation noise E2 


0.006039 0.004196 
0.033567 0.031173 
Balanced estimator 


0.002933 
0.029393 


-0.52 
-0.10 


Poisson noise Epo 
Correlation noise E2 


0.002155 0.001300 
0.009067 0.007152 


0.000803 
0.005812 


-0.72 
-0.33 
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of noise disappears if wc have information on the redshifts 
of galaxies (e.g., though photometric redshifts) and if we 
make use of this data. A second important point is the de- 
pendence of the noise terms on the lens redshift. Figure ^ 
gives a numerical estimate for the various sources of error 
as a function of Zd- AH errors increase with redshift, and 
in particular the error due to source ellipticities is quickly 
increasing. In fact, the lensing signal decreases as 1/{Z), 
thus is strongly dependent on the lens redshift (we recall 
that Z{z) = for z < Zd). Other sources of noise are less 
dependent on the redshift of the lens. It is also interest- 
ing to observe that the largest contribution to Poisson and 
clustering noise comes from the angular dependence [com- 
pare the curve Epo for the unbalanced estimator, which is 
relative to Poisson noise in angular positions and redshifts 
of galaxies, with the curve Spo for the balanced estimator, 
which includes only redshift Poisson noise; similarly for 
the curves E2] . This last point is basically due to the "com- 
pression" of distances in the redshift space. Assuming that 
the redshift is only cosmological, two galaxies at slightly 
different redshifts (say Az — 0.01) can actually be several 
Mpc away, so that only weak correlations are expected 
for such galaxies. On the other hand, their cosmological 
weight is basically the same, and thus no extra noise on the 
shear is expected. We also stress here the different angular 
dependence of the clustering noise, which decreases just as 
l/\6 — 0'\ at large radii (see Sect. |C.4[ ), i.e. it introduces a 
long-range correlation on the map (in contrast, ellipticity 
and Poisson noise correlation decrease very rapidly with 
|0 — as a Gaussian). 

Finally, it is worth to mention a general technique of- 
ten used to estimate errors of weak lensing studies, namely 
bootstrapping. This method is based on the generation of 
many datasets from the observed data to be used for the 
determination of errors. Suppose that we observe N galax- 
ies with ellipticities {e^^"^} and positions j^*^")}. We can 
then construct a new galaxy catalog by randomly draw- 
ing N galaxies from the original catalog, with replace- 
ment. Because of the replacements, we are allowed to use a 
galaxy, which is represented here by {0, e), more that once. 
We can then use the new catalog to perform a weak lens- 
ing mass reconstruction. By repeating the whole process 
several times, we eventually obtain the noise properties of 
the mass map. Bootstrapping is a very simple and rea- 
sonably robust method to obtain error estimates, but un- 
fortunately relies on the assumption of independent data. 
In our specific case, because of the spatial correlation of 
galaxies, we are not really dealing with independent data, 
and thus, at least in principle, the use bootstrapping is not 
fully justified. Moreover, the use of analytical expressions 
provided here for the noise estimate has the advantage of 
showing clearly the different noise contributions and of al- 
lowing us to predict the expected errors without using any 
specific dataset. 



8. Sources with known redshifts 

When some kind of knowledge about the redshift of each 
galaxy is available, for example through photometric red- 
shifts, we can take advantage of the extra information to 
reduce the error on the shear (and thus on the mass map) 
by adopting some improved estimator. 

Suppose that measured redshifts {z*^")} for source 
galaxies are available, together with their expected errors. 
Then, assuming a particular cosmological model, we can 
convert these redshifts into the corresponding cosmolog- 
ical weights If the estimates of the redshifts are 
unbiased and their expected errors small, we can take the 
cosmological weights unbiased, i.e. with some "symmet- 
ric" scatter around the true weights. 

In order to take into account errors on measured red- 
shifts, we introduce pz{z\z), the Baysian probability dis- 
tribution for z. This function gives the probability of mea- 
suring z when the real redshift is z. Note that we assume 
that the probability distribution for z depends only on z, 
the true redshift. Referring to photometric redshifts, we 
observe that, on one hand, this approximation neglects 
some key points, such as the dependence of the error on 
the photometry (and thus luminosity) of the galaxy; on 
the other hand, since, except for the redshift z, we are 
dealing with quantities which can be taken to be uncorre- 
lated with the measured redshift z of a galaxy {9 and e'' 
have clearly nothing to do with z), we can safely simplify 
the discussion using pj(z|z). When performing averages, 
we need to include the new random variable 2. As a result, 
we have to modify Eqs. (^) and ^ with the replacement 

/>oo 

Pz{z) / pz{z\z)dz . (52) 
Jo 

Note that in terms of the probability distribution pz{z\z), 
the condition of unbiased cosmological weights estimates 
can be written as 

(Z) = J pi{z\z)Z{z) dz = Z{z) . (53) 

In order simplify the discussion, in the following we will 
assume that the expected error on the cosmological weight 
is a simple function of the redshift of the galaxy, which we 
write cr^(z): 

1*00 

a^(z)= / p,{z\z)[Z{z)-Z{z)\'dz. (54) 
Jo 

We will often refer to the variance of measured redshifts. 
This quantity is defined as 

/>oo 

= / pz{z)a^{z)dz . (55) 
Jo 

We stress that this variance depends on both the errors 
in redshifts and on the source redshift probability distri- 
bution. To give an example of the precision that can be 
reached by photometric redshifts when accurate photom- 
etry is available, we note that for the Hubble Deep Field 
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North a scatter of the order of 0.07(1 + 2), without outliers 
and bias, has been reached (Benitez 200C). 

We will now investigate two possible uses of the quanti- 
ties {Z^")} for shear estimators. Our starting model will 
be the estimator of Eq. (j^), since the estimator has 
already been proven to have a larger variance. 



adopted in Sect. 4.3, we obtain a vanishing contribution 
to the covariance. In other words, the availability of some 
redshift information completely cancels the noise due to 
the redshift correlation of background galaxies. In deriving 
this result, a key role is played by the assumption that the 
measured cosmological weights are unbiased. 



8.1. Unweighted estimator 

In Sect. I we have seen that the use of a balanced estima- 
tor on the source positions completely eliminates Poisson 
noise due to the galaxy positions and significantly reduces 
the correlation noise. Since estimates for the cosmological 
weights are now available, we can now think of replacing 
the average (Z) in the denominator of Eq. (^) with the 
measured values for the cosmological weight. This leads 
to the new shear estimator 



(56) 



This estimator is called "unweighed" since we are not 
taking advantage of the redshift information to weight 
galax ies (a weighted estimator will be discussed below in 
Sect. B^). Using a technique very similar to the one dis- 
cussed above, we can study the statistical properties of 
this estimator, and in particular its mean and variance. 

Regarding the mean, in the linear approximation this 
estimator behaves exactly as the estimators considered so 
far, i.e. Eq. ( [Tsl ) holds. In the slowly varying shear approx- 
imation, moreover, we have 0{0)) = 7(0). We stress that 
a key point in deriving these results is Eq. (|5^), i.e. the 
condition of unbiased cosmological weights. 

Using the slowly varying shear approximations and in 
case of a vanishing two-point galaxy correlation function, 
the covariance of 7 is 



Covy(7;0,0') 



1 



X exp 



2 

\e- 



1/12 



{Zf 



1-47.(0)7^(0' 
7.:..'(0)7w'(0) 



where tr^ is defined in Eq, 



IS very 
2. In 



cr 



^). Note that Eq. (| 
similar to Eq. (^), with the replacement cr| 
general cr^ is smaller than (t|, since a redshift estimate 
will hopefully provide a better constraint on the galaxy 
cosmological weights than using a straight average for all 
galaxies. We stress, moreover, that cr^ is mainly set by 
the redshift inaccuracies, since it represents the scatter of 
a measured cosmological weight with respect to the true 
value, while depends only on the galaxy redshift dis- 
tribution. Clearly, if no errors on the redshift estimation 
are present we have a| = 0. 

Calculations for the case of a non-vanishing galaxy 
two-point correlation function, not reported here, are 
rather surprising. In fact, using the approximations 



8.2. Weighted estimator 

The use of an unweighted estimator is rather convenient 
for the discussion of the noise properties of the shear, but 
is not optimal. In fact, we can take further advantage of 
the knowledge of the redshifts of galaxies by suppressing 
galaxies which are (or are taken to be) in front of the clus- 
ter and enhancing distant galaxies, which are well affected 
by the lens. A simple investigation shows that the weight 
to use for each galaxy is proportional to Z^"^ As a result, 
we can think to use the shear estimator (see also 
|fc Bertin 1999bD 



Lombard 



7(0) = 



En=ie'"^^^"^W/(0-0(")) 
ELi(^("^)V(0-0(")) 



(58) 



Unfortunately, this estimator turns out to be biased if the 
estimated cosmological weights Z^"^ are not identical to 
the true ones Z'"'. In fact we have 



(7(0)) = 



{Z-) 



Jn 



(59) 



In other words, we will underestimate the shear by a factor 
l+o-|/(Z)2. This, hopefully, is a number very close to one, 
so we will for the moment ignore this effect (for example, 
assuming cr^ < cr^ , this factor differs from 1 by an amount 
smaller than 0.02 for a lens at redshift Zd — 0.3). 

The expression for the covariance of 7 is in general 
rather complicated because of the presence of complex 
combinations of averages of Z and Z. For this reason, 
we do not explicitly write here the general expression. In 
the special case of no errors (so that Z = Z), however, we 
obtain the simple result 



(57) Cov,,(7;0,0') 



1 



An{Z^)pa^^ 



{Z^)a 



w X 



i'ji 



X exp 



2 

|0_0]2 



Vi'^r 7«,»' (0)7j J' (0) 



(60) 



This equation should be compared with Eq. ( |57| ) with 
(t| = 0. We see that the main difference is the factor {Z"^) 
which replaces {Z)"^ . We can thus conclude that, in case of 
complete redshift knowledge, the estimator ( p8| ) performs 
better than the one of Eq. (|5q). Actually, the gain is rather 
low (about 10% for the lens at redshift Zd = 0.3), and one 
could also prefer to keep the unweighted estimator (56) 
rather than the weighted one as long as the lens redshift 
is much smaller than the median redshift of galaxies. 
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8.3. Discussion 

As shown above, one can well take advantage of redshift 
information to reduce the noise of shear estimates. Since 
weak lensing analysis normally deals with thousands of 
galaxies, redshift information is generally available only 
though photometric redshifts, and thus errors on redshifts 
must be expected. 

In general, galaxy redshifts, even with errors, are very 
valuable data: They can be used to severely reduce the 
noise of the mass reconstruction, especially the contribu- 
tion due to galaxy correlations; moreover, they can also be 
used to select background galaxies or, more generally, to 
weight galaxies depending on the lensing signal that they 
carry. On the other hand, some care should be used when 
photometric redshifts are used to select galaxies. In fact, 
errors on photometric redshifts are extremely sensitive to 
the accuracy of the photometry and on the bands used. 
Large errors, clearly, imply some kind of "contamination" 
between foreground and background galaxies and, ulti- 
mately, introduce a bias on the shear estimate [see above 
Eq. g)]. 

In general, we find that photometric redshifts can al- 
most always be used together with the estimator ( ^6|) 
(which, we recall, can be significantly less noisy than es- 
timators not accounting for galaxy redshifts). We think, 
instead, that the estimator (p^), should be used only if er- 
rors on the redshifts are statistically well known and not 
too large or, alternatively, if the lens is at high redshift. 
In fact, for a high redshift lens, we expect a significant 
fraction of foreground galaxies: The use of a "filter," pro- 
vided by the weight Z^") in the numerator of Eq. (|5^), can 
significantly reduce the noise of the shear and thus of the 
mass map. 

9. Conclusions 

In this paper we have presented a careful statistical analy- 
sis of the shear estimators most widely used in weak lens- 
ing mass reconstructions. In particular, we have obtained 
analytical expressions for the bias and the covariance of 
the estimators, and identified the various contributions to 
the noise, as summarized in the following items. 

— Intrinsic spread of ellipticities. This term, which is 
common to all estimators, is basically the only one 



considered by previous investigations (see Kaiser et al, 
1995| ; [Lombardi fc Bcrtin 1998a| ; [van Waerbeke 2000| ) 



Poisson noise. Poisson noise is due to the spread of 
galaxies in positions and redshift. We have shown that 
a simple, balanced estimator is able to eliminate the 
angular component of Poisson noise, thus reducing this 
source of error by a factor ~ 4. If redshift information 
is available, redshift Poisson noise can also be reduced 
using suitable estimators. 

Correlation noise. The clustering of galaxies introduces 
a further source of noise. We have shown that this is 
generally the leading noise contribution, even for the 
balanced estimator. This source of noise, in contrast to 



other sources, does not formally decrease with the den- 
sity of background galaxies. This error can be severely 
depressed if redshift information is available. 

We have also discussed the dependence of these noise 
terms on several relevant parameters, such as the lens red- 
shift, the lens strength, the density of background galaxies, 
and the smoothing length used. Finally, we have tested 
our analytical framework by performing numerical inte- 
grations in a specific model, and by comparing the results 
obtained with our analytical predictions. We have then 
proved that our finite-form expressions for the shear noise 
are correct to within 10%; moreover, we have shown that 
these expressions can be safely taken also as good approx- 
imations for the error on the mass map. 

Acknowledgements. We would like to thank Giuseppe Bertin 
for helpful discussions. 

Appendix A: Linear approximation 

For the study of the properties of the estimator ^ we have 



used a stan dard technique often applied in statistics (Eadie 
^t al. 1971 , see, e.g.). Suppose that X is a multidimensional 
random variable and Y — f{X), with / a smooth function. If X 
has small covariance, we can expand / to first-order and obtain 
this way a simple expression for the mean and the covariance 
of F: 



(Y) ^ f{{X)) , 



Cov,,(F)~C,,,C,y CoY,,,,{X) 



where 



Cij — 



dfi 
dXi 



(A.l) 
(A.2) 



(A.3) 



X^{X) 



The goodness of the first-order approximation can be verified 
by using a second-order expansion. Galculations are particu- 
larly simple for the average value: 



(Y) ^ f{{X)) + - 



■Gov„(X) 



(A.4) 



2 dXidXj 

The second term represents the second-order correction and 
can be used as an estimate for the error made. 



Appendix B: Detailed statistical analysis 

In this appendix we will carry out in some detail a statistical 
analysis for the shear estimator (^ , thus recovering the results 
stated in Sect. 0. 



B.l. Average value 

The average value of the estimator is easily evaluated using the 
linear approximation. In particular, we have 



{N}=p{z} [ wie-0')-/(e')d^ 

Jn 

(D) = p{Z) . 



As a result, using Eq. (A.l), we find 

(7(6*)) = / w{e-e')j{e')(fe' . 

Jn 



(B.l) 
(B.2) 

(B.3) 
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B.2. Covariance 

Calculations for the covariances are more difficult. The scatter 
of the numerator from the average value can be written as 



A 



W'(6>-0)7,(</.)d"<^ 



n 



For the denominator we have 

N 



D~{D) = (Z) J2 



n = l 



w{e-e 



A 



(B.4) 



(B.5) 



Using Eq. (B.4), we can write the covariance Cov (NiNj) as 

Cov{NiN'j-e,e') = 



E 



E 



A 



e^f™' + 7, (e^^M Z^™M (6>' - 6> 



l¥(6»- 0)7,(0) d"0 



{m)\ ;7(™) \ II/ l-fl' _ />(m)\ 



^ |^W(0'-0')7. (</>') d'<^' 



(B.6) 



Similarly to Eq. (^l|), it is convenient here to write 
CoY{NiN!j) = Ki + if 2, where Ki includes terms for which 
n — m, and K2 terms for which n ^ m. 

If n = m, Eq. (^ can be used, leading to 

Ki = pa^Jij [ W{e - 0) W(6»' - 4>) d^0 

J SI 

+ {Z^)p[ 7.(0)7. (</')M^(0-</')W^(0'-</')d'0, (B.7) 

where terms which vanish in the limit A ^ <x have been dis- 
carted. If n 7^ m we have 



f dVW/(6>- 0)7,(0) / dVVK(6>'-0')7j(0') 
n in 

rp.{z)Z{z) Az rp,(z')Z{z')a4> - </>', z') dz' . 
Jo Jo 



(B.8) 

We recall [see Eq. (^^] that this integral involving 
the two-point correlation function has been denoted as 
E[WjiZW'j'jZ'] . Finally we find for Cov(7ViiVj) 

Cov(7V,7V;; 6», 6»') = pa^S.j [ WiG - 0)W(6»' - 0) d^ 



+ {z^)p f 7.(</')7.(0)w^(0 - 4>)w{e' - 0) d20 

+ p^H[Vy7iW7jZ'] . (B.9) 

Covariances for the other terms can be calculated in a sim- 
ilar manner. The final result obtained for Cov{NiD') is 



Cov(iv,L»'; e, 6') = {z)^p [ w{e - cj))w{e' - cj))-f,{4>) d" 

in 



+ {Z)p^E[W'y,ZW'] 



(B.IO) 



The result for Cov{DN'j) is analogous and will not be written 
explicitely. Finally, the last covariance needed is 



Cov{DD';e, e') = {zfp / w{e - ct>)w{e' 

Jsi 

\2 



+ {zy p^=.[ww'] 



(B.U) 



In principle, we could now evaluate the covariance of 7 
using the results obtained so far. In practice, the resulting ex- 
pression is by far too complex to be of any practical use at 
this stage. We thus avoid to write this expression in the gen- 
eral case, but rather we will consider the slowly varying shear 
approximation in the next appendix. 

Appendix C: Slowly varying shear 

We will derive here the slowly varying shear expressions for the 
results obtained in Appendix Again, only the estimator (^ 
will be considered; calculations for the estimator (^) are very 
similar. 



C.l. Average value 

To first order, the avera ge v alue of 7 is just the ratio between 
(iV) and (D) [see Eq. (At|)]. Let us focus on (TV), which is 
given by Eq. (B.l). Inserting Eq. (M) in Eq. (B.l) and using 



the symmetry of W we get 

(iV) = p{z)^{e) [ w{e - e') A'e' = p{z)^{e) . (c.i) 

Jn 



Finally, using Eq. (4..1) we obtain 
(7(0)) = 7(0) • 

C.2. Covariance = 0) 



(C.2) 



Calculations for the covariance can be similarly carried out 
using an expansion for 7 around the mid-point 6 = {6 + 0')/2 
(see Fig. |l|). This is convenient because with this choice several 
linear terms containing partial derivatives of 7 vanish. Hence 
we have 



Covy(7;6»,6» 



7» (0)7.(0) 



(Z)2 ' ' (Z)2 

- 7.(0)7.(0) - 7.(0)7.(0') + 7.(0)7.(0') 
X (W* W)(0-0') 

^7.,.'(0)7.,.'(0) 

X / M/(0-0)W^(0'-0)(0,, -0,0 -^.')d'0. 



(C.3) 

We can apply once more the slowly varying shear approxima- 
tion to this result. In the case where is a Gaussian the final 
result can be written in a simpler form: 



Cov,,(7;0,0') = 



+ 



X exp 



[Z^)cTn- 



47v{Zypal, 
{Zf 



2 

|0-0'|^ 



fi'^j' )7i,i'(0)7.,j'(0) 



(C.4) 
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Fig. C.l. Similarly to Fig. this graph shows that linear 
terms of the expansion of 7 around 6 do not contribute to 
the term C2. In fact, the integrand evaluated at the up- 
per configuration (^1, <p'i) is the opposite of the integrand 
evaluated at lower configuration (^2, 02)- 



i.e. Eq. 



C.3. Check of the linear approximation 

We can explicitely verify the goodness of the first-ord er a 
proximation by checking the second-order terms of Eq. (|A 
These terms are both of the form 



7(0) 



wie')w(e')d'^e' 



7(g) 
pAw 



(C.5) 



The two terms have opposite sign. In other words, the average 
value of 7 obtained above is correct to second order [actually 
it is easy to show that the estimator ^ really is unbiased 
to all orders]. We also note that the relative contribution of 
this terms is of the order l/{pAw) ~ 1/Nw, where, we recall, 
Nw the expected number of galaxies inside the effective area 
of W. Note that 1/Nw is also the expected relative variance 
of the denominator D of Eq. (^). We have thus recovered the 
condition needed to perform the linear approximation, already 
discussed at the end of Sect. 

C.4. Covariance (£, 0) 

Given the difficulties encountered to study the contribution 
to the covariance from the two-point correlation function, we 
will carry out the calculations only in the case where is a 
Gaussian of the form 

Let us focus on a specific term involving the two- point cor- 
relation function, for example the leist term of Eq. (B.9), i.e. 



L^ij = S^W^fiZW'^f'jZ'] [in this section, for simplicity, we will 
use the symbol s Hi, , Hi, an d H to denote, respectively, the last 



z' is very close to z (this property will be repeatedly used be- 
low). On the other hand, all other functions involving redshifts, 
namely Da, Dp (but also pz and Z), are slowly varying with 
z. We thus can approximate ^ as 



^{S,z,z') 



dl{l + zy 



Dl{z)\S\^ + [D'p{z)Yiz'-zr 



(C.7) 



Moreover, since ^ is basically zero for z' 7^ 2, we can take as 
constant the term pz{z')Z{z') in Eq. (p3[), and perform directly 
the integration over z' . We get then 



(C.8) 



d'0Vy(0-0)7,(</)) / dVW(f'-</'')7j(0') 
Jn 

r oc r oc 

/ [p..(z)Z(z)\' dz / ^{4>^4>',z,z')Az' 
Jo Jo 



Using Eq. (C.7) the last integration is trivial: 
^ ' Da{z)D'p{z)\S 



l^arctan 



with 



to = 



D'p{z)z 
'Da{z)\5\ 



(C.9) 



(C.IO) 



In normal conditions to has always a very large negative value 
(unless 2 is exceedingly small; these cases, however, are of no 
interest in weak lensing studies). We thus can safely replace to 
with —00 in Eq. (|g.9|), thus obtaining 



7rd?,(l -I- zy 



^{S,z,z')dz' _ f ^rl' t Mxt 
Da(z)D'p{z)\S\ 



as, 



(C.ll) 



We now can perform the two integrals on G and G' . Since the 
results obtained so far depends only on S — 4>' — (j), we have to 
evaluate two "quasi-convolutions," the "quasi" being due to the 
fact that the integrand also contain terms like 7j(</>) and 7-, (</)'). 
Hencewe can apply the slowly varying field approximation (see 
Fig. C.l), obtaining the integral 



poo p 

l.{G)j,{G) / [p4z)]^[Z{z)]'dz / W^(0-</.)dV 
Jo Jn 

W{G' - 4>')C{(I)' - cl),z)d'^(f> . (C.12) 



Using the properties of convolutions we can recast the two 
integrals on the angular variables in the form 



z)= f ((</>, z)(W *W)(>fi- <P) dV 
Jn 



(C.13) 



where = G' — 6. Using the Gaussian weight function (Q) and 
the result obtained above for C,, we can perform the integration 
in polar coordinates. The result is 



terms of Eqs. ( [B.9[ ), (gTO|), and ( |B.ll[ )]. We will use the ex- 2 
pression (^2|) for the function ^ with the approximate exponent v{if>, 2) = ^ d^jl + z) f 
?7 = 2. We then write ' 2Da{z)D'p{z)gw V' 



exp 



C(^,z, 2') 



dg(l + 2)" 



(C.14) 



D\{z)W+[Dp{z)^Dp{z')y 



(C.6) 



where Jo is the modified Bessel function of first kind. The final 
result for the covariance term is thus 



where 5 = — 0' must be expressed in radians. In this equa- 
tion Da{z) and Dp(z) are, respectively, the angular diameter 
distance and the proper distance of an object at redshift 2. 
The length-scale of ^ is set by do, and typically is a few Mpc. 
As a result, the two-point correlation function vanishes unless 



2gv 



{i + zr[p^{z)z{z)\' 

Da{z)D'p(z) 



dz 



(C.15) 
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0.2 0.4 0.6 0, 



1.2 1.4 



As a result, we get a much larger error, since the effect of 
cancellation of terms clearly shown by Eq. (C.21) is no longer 
present. In particular, for the typical case considered here, we 
expect a correlation about 19 times larger, i.e. an error about 
4 times larger. 



Appendix D: Calculations for estimators with 
source redshifts 

In this appendix we briefly sketch the calculations leading to 
the results described in Sect, bl 



Fig. C.2. The dimensionless quantities as a function 
of the lens redshift Zd- 



The function rj{x) = 7o(a::) exp(— x) represents the leading an- 
gular dependence. The last integral cannot be recast in a simple 
form and must be evaluated numerically (it also depen ds on 
the cosmological model). We stress that, since S in Eq. (C.6) 
has been taken to be expressed in rad ians, we similarly need 
to express aw in radians in Eq. (C.15). 



Calculations for the other integrals involving the two- 
point correlation function, namely Hi and H, are very similar. 
Defining the dimensionless quantities 



T 



dz 



2{Z}'= Da{z)D'p(z) 
the various two-point correlation terms can be written as 



(C.16) 



aw 

aw V 80-4. 



Ti , 



T2 



aw 



V 



8^w 



To 



(C.17) 
(C.18) 
(C.19) 



In summary, the contribution to the Cov(7) from a non- 
vanishing two-point correlation function is 



C2 = — [t27.(0)7,(0) -f To7.(0)7,(0') - Ti7,(0)7,(0) 
aw L 



Ti7.W7.(0'; 



\0-0'\ 



(C.20) 



Note that ri{x) oc l/y^ for large x, and thus this contribution 
to the covariance is the leading term at large separations \(f)\. 

In the case — 0' this expression can be written in the 
more compact form 



/■ 2 + iO — ^ll .„^ 

aw 



(C.21) 



For a lens at redshift Zd = 0.3 we have T2 ~ 2.8267 • 10"^, 
Ti ~ 2.6513 ■ 10"^ and Tq ~ 2.6246 • 10"^ so that T2 + Tq - 
2Ti ~ 1.488 ■ 10"^. A plot of the qua ntiti es Tfc as a function 
of the lens redshift is provided in Fig. C.2. 

For Eq. (M) we obtain similar results, but now only the term 



ij contributes to the covariance. More specifically, we find 



C2 = ^7.(0)7.(0)^ 

aw 



>W 



(C.22) 



D.l. Unweighted estimator 

Let us call A'' and D the numerator and the denominator of the 
r.h.s. of Eq. (p6|), respectively. Assuming unbiased measured 
cosmological weights, we easily find 



{N)^p{Z) ( W{0-0')^{0')A^e\ 
Jn 

(D) = P{Z) , 



(D.l) 
(D.2) 



so that Eq. (^6|) is recovered. 

Regarding the covariance, we find after some calculations 

Cov{N,D'-0,0') = 

{Z'')p f W(6»-</))H/(6>'-0)7.(</))d'0 
Jn 

+ p'-lW-y^ZW'Z'] , (D.3) 

Cov(_DD';0,6>') = / W{0-(i))W{0' -(t>)A^(l) 

Jn 

+ p'^E[W ZW' Z'] . (D.4) 



Note that, since the num erato r is formally identical to the nu- 
merator of Eq. Eq. (B.!;) still holds. Differences with re- 
spect to t he analogous terms calculated above in Eqs. ( |B.10| ) 
and (B.ll) are the presence of (Z^) and {Z^) factors instead of 
{Z)'^ and the different form of the correlation terms, which now 
include always the redshift Z. Substituting these expressions 
in Eq. ( ^ ) and using the slow varying shear approximation we 
obtain Eq. (||). 



D.2. Weighted estimator 

The slowly varying shear approximation applied to the numer- 
ator N and denominator D of Eq. (|5q) gives 



{N) = p{z) [ w{0-0'hi0')d^e' , 

Jn 

(D) = p{Z) . 

These equations lead immediately to Eq. (|5£|). 



(D.5) 
(D.6) 
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Regarding the covariances, we find 

CoY{N,N--e,6') = al{z'^)5^, ( w{e -e")w{e' -e")A^e" 

Jci 

+ {z'^z^)p ( 'y^{e")'yJ{e")w{e^e")w(e' ~e")d^e" 
Jit 

+ p'~[W"i:z"Z"W"'-i'^'Z"'Z"'] , (D.7) 
Cov{NiD'-9,e') = 

+ / w{e -e")w{e' ~e")j,{e")<fe" 

Jn 

+ p'~[W"-f':z"Z"[Z"'fw"'] , (D.8) 

Cov{DD'- e, e') = {z^)p I wie - e")w{e' - e") A^e" 

Jn 

+ p^E[W"iZ"fW"'{Z"'f] . (D.9) 
These expressions can be used to derive Eq. (|60|). 
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